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ABSTRACT 

Observations in the sub-THz range of large solar flares have revealed a mysterious spectral 
component increasing with frequency and hence distinct from the microwave component com- 
monly accepted to be produced by gyrosynchrotron (GS) emission from accelerated electrons. 
Evidently, having a distinct sub-THz component requires either a distinct emission mechanism 
(compared to the GS one), or different properties of electrons and location, or both. We find, 
however, that the list of possible emission mechanisms is incomplete. This Letter proposes a 
more complete list of emission mechanisms, capable of producing a sub-THz component, both 
well-known and new in this context and calculates a representative set of their spectra produced 
by a) free-free emission, b) gyrosynchrotron emission, c) synchrotron emission from relativistic 
positrons/electrons, d) diffusive radiation, and e) Cherenkov emission. We discuss the possi- 
ble role of the mechanisms in forming the sub-THz emission and emphasize their diagnostics 
potential for flares. 

Subject headings: Sun: flares — acceleration of particles — turbulence — diffusion — Sun: magnetic 
fields — Sun: radio radiation 



1. INTRODUCTION 

Solar flares, being manifestations of prompt 
energy releases, produce electromagnetic radiation 
throughout the entire spectrum of electromagnetic 
emission from rad i o wa v es to gamma rays (e.g . 
iDennis fe Schwartz! 119891 : iBrown fc Kontari 120051 ). 
Although the radio, UV, X-ray, and gamma-ray 
emission is generally well observed with high 
resolution from huge number of events, the sub- 
THz radiation has only recently been observed 
from a few large events, at a sm a ll number of 



frequ enc ies dKaufmann et al.l l2001 | : | Trottet et al 
20021: iLiithi et all I2004al lbt iKaufmann et al 



2001 iCristiani et al.l 1 2001: iTrottet et all 120*01 



Kaufmann et all l2009blfa l: iSilva et all l2007f ). The 



available observational tools are unable to measure 
polarization and are clearly insufficient to provide 
detailed spectral and positional information about 
the sub-THz bursts. Nevertheless, available sub-THz 
observations have already provided us with puzzling 
questions which have yet to be answered. 

The observations suggest that, on top of quiet Sun 
emission, at least two kinds of sub-THz emission can 
be produced. The first kind looks like a natural exten- 
sion of the microwave spectrum at higher frequencies 
and so can reasonably be interpreted as synchrotron 
radiation from accelerated electrons, which are also 
responsible for microwave and hard X-ray emission. 
The second kind looks like a distinct spectral com- 
ponent rising with frequency in the sub-THz range 
in contrast to the microwave spectrum, which falls 
with frequency. The origin of this component is un- 
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clear; there is no consensus about the possible emis- 
sion mechanism producing it, moreover, it is ambigu- 
ous whether the emission is of thermal or nonthermal 
origin. 

The main observational characteristics of this com- 
ponent arc: relati vely large radiation peak flux of the 
order of 10 4 sfu (|Kaufmann et all 120041 ): radiation 
spectrum rising with frequency F(f)(xf s ; spectral in- 
dex varying with time within 5 ~ 1 ... 6; sub-THz 
component can display a sub-seco nd time variabil- 
ity wit h the modulation about 5% (jKaufmann et all 
l2009aTi : the source size is believed to be less than 
20", however, this conclusion is based on a multi- 
beam observation of a few antennas with ~ 4' reso- 
lution each, rather than on true imaging; therefore, 
the size estimate must be considered with c aution. In 
addition, observations (jLiithi et al.l l2004aT TbT) suggest 
the existence of both compact ~ 10" and extended 
~ 60" components, and source sizes increasing with 
frequency. 

The emission mechanisms proposed so far to ac- 
count for the sub-THz component are thermal free- 
free emission, GS emission from flare accelerated elec- 
trons, and synchrotron emission from n uclear decay- 
gener ated relativistic positrons (e.g. iNindos et al.l 
I2008L as a review). None of these mechanisms 
is readily consistent with the full list of the sub- 
THz component properties and/or with available con- 
text observations at other wavelength. Below we 
analyze these options and also consider two other 
emission mechanisms capable of producing a spec- 
trum rising with frequency — diffusive radiation in 
Langm uir (DRL) waves (e.g. iFleishman fc Toptvginl 
200 7313) and Vayilov-C herenkov emission (e.g., 



Bazvlev fc Zhcvago ll987|) from chromospheric layers, 



which make the list of options more complete. Calcu- 
lating the spectrum from all of the above mentioned 
models, we establish the range of main source param- 
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eters for each model emphasizing the strengths and 
weaknesses of each emission mechanism. 

2. FREE-FREE EMISSION 

Perhaps the simplest example of a radio spec- 
trum rising with frequency is optically thick ther- 
mal free-free emission. Having a rising spectrum 
from a compact (< 20") source requires obviously 
that the source is relatively dense (n e >10 n cm~ 3 ) 
and hot ( T e >10 MK). From available X-ray observa- 
tions (e.g. iDennis fe Schwartdfl989l ; iKasparova et al.1 
we can exclude the option of a source that is 
simultaneously dense and hot, say n e ~ 10 12 cm~ 3 
and T e ~ 10 MK, EM = n 2 e V - 3 x 10 51 cm" 3 , 
since a plasma with the corresponding emission mea- 
sure and temperature would produce stronger X- 
ray e mission than is actually observed from solar 
flares (|Dennis fe Schwartz! [19891 : lEmslie et all [20031 : 
iBrown et all 120071) . We cannot exclude, however, 
dense plasmas of lower temperature, T e ~ 1 MK; such 
plasmas do not contradict to the current X-ray and 
UV observations. 

Figure [T] shows free-free spectra from overdense 
sources with temperature around 1 MK. Evidently, 
having a flux density above 1000 sfu level requires 
thermal electron number density above 10 12 cm~ 3 
or/and the linear size of the source above 20". Al- 
though available size estimates do not favor sources 
larger than 20", a firm conclusion about that must 
await true imaging observations at the sub-THz 
range. Such large emission sources, tens of arcsec- 
onds, are consisten t with at least some observations 
(ILuthi et al.l[2004aH bT>. 

It should be noted that the spectral index of 
a spatially uniform source of free-free radiation 
cannot be larger than 2, while the observations 
suggest larger values, especially during the ini- 
tial ph ase of the sub-THz bursts (jKaufmann et all 
2009b). These larger values can be reconciled 
with the free-free mechanism if we allow for a 
cooler dense absorbing layer in the line of sight 
between the source and observer, providing an at- 
tenuation factor exp(— T//), where T//ocnro e (l — 
exp(-/i//fcT))T- 1 / 2 /- 3 oc mi e T- 3 / 2 /- 2 is the free- 
free optical depth; such a non-uniform source is capa- 
ble of producing frequency spectra with a very sharp 
low frequency edge compatible with observations. 

Finally, wc emphasize that the free-free emission 
mechanism can display temporal variability. For 
sausage mode loop oscillations, for example, B = 
Bq(1 + m cos(2wt / P)) , with a period P and modu- 
lation amplitude m < 1, we can calculate all relevant 
parameter variations and, thus, make firm predic- 
tion about the free-free emission modulation ampli- 
tude in the optically thick and thin regimes: F t hick °c 
(1 + f cos(27rf/P)) and F thin oc (1 + C os(27rf/P)). 
Therefore, the optically thick and thin part of the 
free-free emission oscillate in phase. However, the 
modulation amplitude in the thin regime is four times 
larger than in the thick regime. 

We can conclude that the free-free model can pro- 
vide sufficient flux and spectral index consistent with 
observations of the sub-THz component, although the 



source density must be somewhat high compared with 
standard coronal loop densities. 

3. GYROSYNCHROTRON EMISSION FROM A 
COMPACT SOURCE 

Flare generated electrons gyrating in a magnetic 
field produce strong gyrosynchro tron emission across 
a broad range of freq uencies (e .g. Basti an et al.l fl998: 
iNindos et aHl2008| ). iSilva et all (|2007j) have demon- 
strated that a rising optically thick GS spectrum in 
the sub-THz range requires a large magnetic field 
in the emission source; typically, this implies also 
a large plasma density. Therefore, we include the 
free-free absorption and emission along with the GS 
contribution in our calculations. We have investi- 
gated the corresponding paramet er space and con - 
fidently confirm the conclusion of ISilva et al.l (|2007f ) 
that only a very compact source ~ 1" with a very 
strong magnetic field B > 2000 G, in which all 
flare accelerated electrons (N e [> 50 keV] = 5 • 10 35 , 
n e [> 50 keV] ~ 10 12 cm~ 3 ) simultaneously reside, 
could be consistent with the spectrum of sub-THz 
flare component, Figure [2] The number of required 
electrons is typical for a lar ge (GOES X-class) solar 
flare (e.g. IBrown et alJl2007f ). 

Although this model is consistent with a small 
source size, it requires more than extreme source pa- 
rameters. In addition, this model cannot account for 
large spectral indices (> 3) of the sub-THz compo- 
nent. The Razin effect in a dense plasma with a 
magnetic field around 800 G was propos ed to form 
a stee per spectrum in the sub-THz range (jSilva et al.1 
l2007f ). However, as seen from Figure [U the flux den- 
sity for B = 800 G is much lower than observed (in 
fact, the free- free absorption additionally reduces the 
flux here). A larger magnetic field (to increase the 
flux level) will require a proportionally larger ther- 
mal electron density to keep the Razin effect in place, 
which will further increase free-free absorption in the 
source, so the observed flux level will not be met un- 
der conditions of strong Razin effect, Figure [2f b) . An 
increase of the source size above 2" with the same 
total number of fast electrons, magnetic field and 
thermal electron density results in a spectrum totally 
dominated by the free-free contribution — the model 
considered in the previous section. 

We note that GS emission model can easily ac- 
count for the observed temporal pulsations of the sub- 
THz component via, e.g., loop oscillations at large 
time scales and/or fluctuat i ons of fast electron in] ec - 
tion (|Kiplinger et all 119841 : lAschwanden et al.lll998l ) 
at various time scales. In addition, good temporal 
correlation with high energy hard X-rays and gamma- 
ray continuum will be rather natural. Nevertheless, 
we feel that unrealistically extreme GS source param- 
eters are needed to reconcile the model with observa- 
tions of sub-THz emission from large flares. In con- 
trast, weaker sub-T Hz fluxes Ff ~ 100 s fu observed 
from M-class flares (jCristiani et al.ll2008f ) may agree 
with the GS emission mechanism. 

4. SYNCHROTRON EMISSION FROM 
RELATIVISTIC POSITRONS 
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Fig. 2. — (a)Radio spectra produced by GS plus free-free contributions from a uniform source with a size of 1" for n e = 8-10 12 cm -3 
and B = 800 . . . 4400 G. (b) Razin-suppressed GS spectra with the Razin frequency 200 GHz plus the free-free component. 



Flarc-accclcratcd ions of tens McV or above col- 
liding with thermal ions in dense layers of the solar 
atmosphere trigger nuclear reactions, with relativistic 
positr ons being one of the product s of such interac- 
tions (jLingenfelter fe Ramatvl Il967l ; iKozlovskv et al.1 
2002]). These relativistic positrons are capable of pro- 
ducing sub-THz synchrotron radiation with a peak 
frequency around /peafc~/se7 2 , where fse is the elec- 
tron gyro frequency and 7 is the Lorenz-f actor, as sug- 
gested by iLingenfelter fc Ramatvl (|1967f ). 

Although the synchrotron peak frequ ency can eas- 
ily fa ll into the sub-THz or THz range (jTrottet et al.l 
I2008T) considering typical magnetic field (B ~ 1000 G) 
and Lorenz-factors 7 ~ 20, the low-frequency syn- 
chrotron spectrum, Fjoc/ 1 / 3 , Figure a) , is incon- 
sistent with the whole range of the observed spectral 
index values unless we adopt an additional absorp- 
tion mechanism leading to a sharper spectral shape. 
The only viable absorption mechanism is free-free 
absorption, which again requires high plasma den- 
sity, ~ 10 12 cm~ 3 , Figure [UJb). Even though a 
combination of the synchrotron radiation from rel- 
ativistic positrons and free-free absorption is capa- 



ble of producing the required spectral shape, the cor- 
rect flux density requires more relativistic positrons 
than seems to be available from nuclear interactions 
even in large flares. The total number of energetic 
protons above 30 MeV for large solar flares as de- 
duced from RH ESSI observation s is in the range 
10 29 - 10 33 (e.g. IShih et al.l[200l . which is compa- 
rable to the instantaneous number of positrons with 
Lorentz factor, 7 = 20 required to explain observed 
sub-millimeter fluxes. The total thick-target positron 
yield is about 1 0~ 2 — 10~ 4 per proto n of energy 
above 10 MeV (jKozlovskv et all 12004 ). so the to- 
tal number of positrons produced in a flare should 
be < 10 31 . Since the positron lifetime in a dense 
positron-production site is likely to be less than du- 
ration of a flare, the instantaneous positron number 
of ~ 10 31 is difficult to achieve, unless we allow a 
significant fraction of relativistic positrons to escape 
to and then be trapped at a tenuous coronal part of 
the flaring loop, where the positron lifetime is much 
longer. 

Synchrotron radiation from relativistic positrons is 
easily consistent with small sizes if the emission re- 
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Fig. 3. — Radio spectra produced by synchrotron radiation from relativistic positron plus free-free contribution from a uniform 
cubic source with a linear size of 20" for the total instantaneous positron number JV e + = 5 X 10 29 ... 5 X 10 31 , with energy 7 = 20 
(~ 10 MeV), magnetic field B = 1000 G, and the thermal electron density n e = 4 ■ 10 10 cm -3 (a) and n e = 4 • 10 11 cm -3 (b), and 
T e = 1 MK. 



gion is placed in footpoints of a flaring loop (< 7") 
or in a moderate-size (< 20") coronal flaring loop; 
such sizes are consistent respectively with hard X- 
ray m easurements of footpoint sizes (jKontar et al.l 
200cf) and radio measure ments of the coronal loop 
sizes (|Bastian et al.|[l998() . 

5. DIFFUSIVE RADIATION IN LANGMUIR 
WAVES WITH LARGE WAVELENGTHS 

Unlike the emission processes discussed above, 
which require dense plasma to produce a sufficiently 
fast-growing spectrum with frequency, there is a class 
of wave up-scattering processes, for which a high 
plasma density is not essential. Inverse Compton 
scattering of gyro-synchrotron emitted photons is 
easy to estimate to be inefficient in solar flare con- 
ditions because the energy density of electromagnetic 
emission is much lower than that of the magnetic field, 
and r;c<l- Nevertheless, similar processes, in which 
the role of low-frequency photons is played by plasma 
waves can be relevant. 

Indeed, a flare-accelerated or decay-produced rela- 
tivistic charged particle moving through the plasma 
with the turbulent waves experiences random Lorenz 
forces and so diffuses in space. Accordingly, the cor- 
responding ra diative processes a re called diffusive ra- 
diations (e.g., [Fleishman 20061). Generally, the dif- 
fusive radiations do not produce a sp ectrum rising 
with frequency (e.g., iFleishmanl 12006) . An exclu- 
sion is DRL for long-wavelen gth Langmuir waves 
(jFleishman fe Toptvginl 2007a. b). In this case the 
spectrum peaks at frequency 2f pe j 2 , where f pe is the 
plasma frequency, Figure 0J The spectral index can 
be as large as 2. The time variability is expected due 
to nonlinear dynamics of the Langmuir waves, which 
are known to o scillate under nonlinear wave- wave in- 
teractions (e.g. lKontar fc Pecseli|[2002T) . 

Although this emission process is attractive, it is 



unclear whether there are sufficient levels of Lang- 
muir turbulence energy density and relativistic elec- 
tron/positron numbers to produce the required flux 
densities. Langmuir waves are effectively generated 
at shorter wavelengths A = 2irv/ujp e < 2Trc/u> pe , at 
the typical resonance phase vel ocities with the elec- 
tron beam v = (0.2 — 0.6)c (e.g.. lKaplan fc Tsvtovichl 
IT^ : IKontanl200l . The refore, non-resona nt genera- 
tion of the plasma waves (IDieckmannll2005f) or nonlin - 
ear wave-wave processes dKaplan fc TsvtovicIj[l97l 
transferring energy from small to large wavelength are 
also implied for this model. 

Typically, a large level of Langmuir turbulence 
would result in strong plasma radio emission due to 
either coalescence of the plasma waves or scatter- 
ing/decay into transverse electromagnetic waves. In 
the DRL model considered here, however, most of 
the Langmuir turbulence energy resides in long wave- 
length Langmuir waves, whose wave vectors are small, 
k <C LUpe/c. These long wavelength waves cannot co- 
alescence into transverse waves because the conserva- 
tion of momentum cannot be fulfilled in this three- 
wave process. Only a minor part of the energy from 
Langmuir turbulence with k > o; pe /c will be involved 
in the plasma radiation production, resulting in a 
modest level of coherent decimeter emission if any, 
which does not contradict the context radio observa- 
tions. 

To round up the discussion involving the plasma 
waves, we note that all nonlinear plasma processes in- 
cluding Langmuir waves and non-relativistic particles 
are of little interest here since the fundamental and 
harmonic plasma frequencies are well below 30 GHz 
in bo th chromosphere and corona (jAschwanden et alJ 
I2002T ). 

6. VAVILOV-CHERENKOV RADIATION FROM 
CHROMOSPHERIC LAYERS 
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(a) (b) 

Fig. 4. — DRL spectrum produced by relativistic electrons/positrons in long-wave Langmuir turbulence X^2irc/u) ve . (a) Depen- 
dency on Langmuir wave energy density, (b) The spectra for different Lorentz factors. 
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Fig. 5. — (a) A model of plasma dielectric permittivity with molecular line contribution included; (b) Vavilov-Cherenkov radiation 
produced by fast electrons with a power-law distribution over the velocity — blue and green curves; the red curve is for e(ui) = 
1 +u) 2 /u)q, i.e., without standard plasma contribution. 



Vavilov-Cherenkov radiation is produced by any 
charged particle moving faster than the correspond- 
ing speed of light in a medium. In fully ionized 
plasma the high-frequency dielectric permittivity is 
less than unity, e(oj) < 1. Therefore, the phase veloc- 
ity of electromagnetic waves cj ^Je{uj) is larger than 
the speed of light, c, and the Cherenkov emission 
docs not occur. However, the chromosphcric gas is 
only partly ionized; and there are numerous atoms 
and molecules whose quantum transitions can make 
positive contribution to the dielectric permittivity, 
making Cherenkov radiation possible at certain fre- 
quency windows. Charged particle of velocity v emits 
Cherenkov emission in the medium, where its dielec- 



tric permittivity s{uj) is such that v > cj \J e(yi). The 
dielectric permittivity of gases is only slightly more 
than unity in the optical range [e.g., Hydrogen gas 
has (e — 1) ~ 2 • 10~ 4 ], so only highly relativistic par- 
ticles with v > (1 — 10~ 4 )c can emit. The situation is 
however very different in other frequency ranges, near 
transition frequencies of atoms or molecules: e(u) 
goes up and down being considerably larger than one 
at certain frequency windows, which allows even sub- 
relativistic electrons to emit Vavilov-Cherenkov radi- 
ation. 

Each atomic or molecular quantum transition 
makes a contribution to the dielectric permittivity of 
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the form of 5s 



47rn e e 



where 



m e (w£ m -w 2 )+ir„ m w ' 

iS nOT is the oscillator strength of the transition, uJ nn 
is the transition frequency, T nm is the transition de- 
cay constant. The resulting dielectric permittivity 
accounting for plasma (free electrons) and molecular 
contributions is s{lo) = 1 — ^p e /^ 2 + J2 & £ nm- The ex- 
act spectroscopic permittivity depends on the chem- 
ical composition of the chromosphere with the sum 
over all excitation states of corresponding molecules, 
which is beyond our intention. Here we perform order 
of magnitude estimates. We assume that there are 
many atomic/molecular transitions capable of mak- 
ing a contribution to the dielectric properties of gas. 
Then, we note that the Cherenkov radiation spec- 
trum rises with frequency only if the dielectric per- 
mittivity raises with frequency. For example, the en- 
ergy level distribution between the rotational levels 
of chromospheric molecules can be assumed to have 
a Boltzmann distribution with the gas temperature 
T, therefore population densities will have maxima at 
energy levels above the THz range increasing to about 
ksT. Thus, we adopt the mean molecular contribu- 
tion to the dielectric permittivity to have the model 
form w 2 /ol>q, where u>q is an unknown constant. 

Assuming a power-law spectrum of fast electrons 
n e (v) = ANqVq v < c, where Nq is the to- 

tal number of electrons with velocity above minimum 
velocity vq, (3 is the spectral index, and A is a dimen- 
sionless normalization constant of order of unity, the 
Cherenkov emission yields flux at Earth: 



Ff = 



10 19 (2ir) 2 AN Q e 2 fvg 



(3-1 



2 £ (/3-2)/2 l 

f3((3 -2) + Je 







(sfu). (1) 



Figure [5(a) displays the model dielectric permittiv- 
ity rising with frequency as described, while Fig- 
ure [5jb) shows the corresponding Cherenkov spectra. 
The spectrum shape and the flux density level allow 
this emission to be easily reconciled mechanism with 
observations for typical numbers of accelerated elec- 
trons. The model flux density is much larger than 
observed, being 5 x 10 7 sfu at 400 GHz, which allows 
us to relax the number of radiating electrons or/and 
add some free-free absorption to the model. Given 
that the emission is from compact footpoints < 10", 
where energetic particles interact with the chromo- 
sphere, we conclude that the Cherenkov emission is 
fully consistent with the observations of the rising 
sub-THz component of large solar flares. The time 
variability requires corresponding fluctuations of the 
electron distribution function, similar to sub-second 
variations of microwave GS or HXR radiations. 

7. DISCUSSION AND CONCLUSIONS 

We have analyzed a number of emission processes, 
which are capable of producing radiation at the sub- 
millimeter wavelengths: free-free emission, gyrosyn- 
chrotron and synchrotron processes, diffusive radia- 
tion, and Vavilov-Cherenkov process. Having in mind 



the characteristic ranges of parameters of a solar flare 
region, we calculated the spectra for each process dis- 
cussed in the paper. Although current observations 
at two frequencies only do not provide us with de- 
tailed spectral shape, we estimate the range of flux 
levels and spectral indices for each model and com- 
pare them with the observations. 

It is likely that the sub-THz emission originates 
from more than a single source and more than one 
mechanism is involved. Free-free emission is a plausi- 
ble candidate in many cas e s, at lea st for large sources 
observed by iLiithi et all (|2004allbh . The free-free 
emission is clearly always present, so other mecha- 
nisms build additional contributions on top of the 
free-free component. Gyrosynchrotron/synchrotron 
emission is likely to play a role in moderate events 
and also as the extension of normal microwave bursts, 
which falls with frequency. The role of DRL is less 
clear since the level of long-wavelength Langmuir 
waves is yet unknown in flares. Finally, the Vavilov- 
Cherenkov emission from compact sources located at 
the chromospheric level seems to be a plausible pro- 
cess to account for the rising with frequency sub-mm 
component of large flares. Indeed, the dielectric per- 
mittivity of partially ionized plasma is known to fluc- 
tuate around unity due to atomic and molecular tran- 
sition contributions. When the flare-accelerated elec- 
trons are present, they will produce Cherenkov ra- 
diation at all frequency windows (including IR, viz, 
and UV bands) where the dielectric permittivity is 
(even marginally) above unity, giving rise to a ra- 
diation spectrum raising with frequency for (mean) 
dielectric permittivity increasing with frequency. 

The available constraints on sub-THz source sizes, 
time variability, and spectral index are not yet fully 
reliable particularly because the ground-based sub- 
THz observations are extremely difficult to calibrate 
due to strongly variable atmospheric opacity. Fur- 
ther progress in understanding the physics of sub-THz 
emission from flares requires observations a) with a 
more complete spectral coverage at the sub-mm range 
(preferably well-calibrated space-based observations) 
and b) polarization measurements. Nevertheless, the 
sub-THz spectral window can be extremely informa- 
tive, e.g., for diagnostics of the chromospheric chem- 
ical composition if the role of the Vavilov-Cherenkov 
emission is confirmed. This calls for a new project 
mounting a sub-THz receivers/interferometers, com- 
bining good sensitivity with high spatial resolution, 
on a space mission, complementing on-going efforts 
in the microwave range. 
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